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ASSmCT 


Temperature  studies  were  conducted  during  both  simulated  and  human 
occupancy  tests  on  the  NBDL  100  man  shelter.  The  temperature  of  the 
shelter  remained  comfortable  during  l4  days  of  human  occupancy.  The 
ventilating  air  movement  through  the  shelter  designed  for  l6  cfm  per 
occupant  was  adequate  to  reduce  the  carbon  monoxide  concentration  from 
heavy  smoking  to  below  safe  tolerance  limits  and  It  removed  69  ^  of  the 
heat  generated  by  the  occupants  with  only  3I  $  of  the  heat  being  dissi¬ 
pated  through  the  walls.  The  shelter  temperature  varied  between  7^  and 
90^  with  an  outside  air  temperature  of  ^0  to  94°F  during  the  simulated 
occupancy  test  and  between  70  and  82^?  during  the  human  occupancy  tests 
with  an  outside  air  temperature  of  36  to  67°!.  Poor  distribution  of 
the  inlet  air  resulted  In  uncomfortably  cool  areas  In  the  front  of  the 
shelter  during  the  human  occupancy  test. 

A  prellmlnazy  analysis  was  conducted  on  a  thermal  analyzer  to 
detezmlne  Its  potential  as  a  tool  for  predicting  shelter  temperatures 
under  various  climatic  conditions.  This  analysis  showed  only  about 
10  $  deviation  from  the  measured  temperature  reported  here. 
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SUMMARY 


The  Problem 

To  determine  the  adequacy  of  the  NRDL  100  man  shelter  design  to 
maintain  comfortable  living  conditions  for  at  least  l4  days  of  human 
occupancy . 

Findings 


Under  the  conditions  of  the  test  the  temperature  and  humidity  of 
the  shelter  remained  within  comfortable  limits  for  l4  days  of  human 
occupancy  test  except  for  cool  drafts  at  the  front  of  the  shelter  due 
to  poor  distribution  of  the  inlet  air.  The  average  heat  output  by  the 
occupants  was  485  Btu/hr  per  man.  The  ventilating  air  movement  through 
the  shelter  of  l6  cfm  per  occupant  was  adequate  to  reduce  the  carbon 
monoxide  concentration  due  to  heavy  smoking  to  below  safe  tolerance 
limits.  It  removed  69  i>  of  the  heat  generated  by  the  occupants  with 
only  31  i>  of  the  heat  being  dissipated  through  the  walls.  The  effective 
temperature  in  the  shelter  varied  between  74  and  80°  during  the  simu¬ 
lated  occupancy  run  when  the  ambient  air  temperature  varied  between  50 
and  94 °F.  During  the  human  occupancy  test  the  effective  temperature 
was  70  to  82°  with  the  outside  air  temperature  varying  between  36  and 
67OF. 


A  preliminary  study  on  a  thermal-electrical  analyzer  at  Oklahoma 
State  University  showed  approximately  10  ^  variation  from  the  measured 
conditions.  This  Indicated  that  reliable  studies  can  be  made  in  the 
laboratory  for  predicting  the  environmental  conditions  in  this  shelter 
located  In  various  climates. 
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CHAPTER  1 


INTRODUCTION 


The  protection  of  the  nation  from  the  hazardous  radioactive  fall¬ 
out  of  a  nuclear  war  has  been  given  much  consideration  in  recent  years. 
The  most  common  considerations  are  the  protection  of  Individuals  in 
family  shelters  and  the  protection  of  large  groups  in  community-built 
shelters . 

A  prototype  100 -man  shelter  vas  constructed  by  NRDL  in  1959  at 
Camp  Parks  approximately  4o  miles  east  of  San  Francisco.  A  variety  of 
experiments  have  been  perfonned  in  this  shelter  to  establish  adequacy 
of  the  ventilation  in  maintaining  comfortable  temperature  and  humidity, 
to  stxidy  the  psychological  effects  of  confinement  in  a  group  shelter, 
to  study  the  need  for  closure  during  mass  fires,  and  fallout  ingress 
e3q>eriraents  to  determine  the  need  for  filters  in  the  ventilation  system. 
The  results  of  these  studies  are  abstracted  in  reference  (l). 

In  the  design  of  a  community  shelter  it  is  essential  that  the  cost 
be  kept  to  a  minimum,  vhether  these  shelters  are  to  be  financed  by  the 
government  or  by  private  enterprise.  The  problem  of  building  a  shelter 
at  lowest  possible  cost  automatically  requires  that  both  the  floor  area 
and  the  gross  volume  per  man  be  reduced  to  a  minimxun.  To  maintain 
physical  comfort  in  a  shelter  under  these  conditions,  sufficient  venti¬ 
lating  air  must  be  provided  to  carry  off  the  heat  and  moisture  generated 
by  the  occupants.  These  factors  were  weighed  against  cost  in  the  design 
of  the  U.  S.  Naval  Radiological  Defense  laboratory  100  man  shelter  and  a 
blower  was  selected  which  was  expected  to  maintain  a  comfortable  envir¬ 
onment  at  minimum  cost.^ 

The  methods  available  for  theoreticeO.  analysis  of  the  heat-transfer 
conditions  in  the  shelter  were  limited  so  it  became  necessary  to  con¬ 
duct  temperature  studies  during  actual  human  occupancy  tests.  As  a 
precautionary  measure  prior  to  conducting  these  tests  a  study  was  made 
with  100  devices  which  simulated  human  occupants. 
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CHAPTER  2 


EXPERIMEmL  CONDITIONS 


2,1  SHEI/EER  COHSTRUCTEON 

Figure  1  Is  an  artist's  sketch  shewing  the  basic  design  features. 

It  Is  a  galvanized  multi -plate  steel  arch  type  structure  (Standard  Navy 
Ammunition  Hut)  with  a  floor  area  of  25  ft  x  lj-8  ft.  It  is  burled  with 
the  top  flush  with  the  ground  level  and  then  covered  with  a  3  ft  mound 
of  soil.  It  has  a  capacity  of  100  people  with  a  floor  space  of  12  sq. 
ft  per  man  and  a  gross  volume  of  117  cu  ft  per  man.  It  has  tiers  of 
bunks  on  each  side  with  a  sleeping  capacity  of  These  bunks  can  be 
rolled  up  to  provide  more  space  in  the  daytime.  Air  intake  openings 
are  located  In  the  entzy-way  at  the  top  of  the  stairs  as  shown  In  Fig. 

2.  Fresh  air  is  drawn  in  with  a  blower  located  outside  the  Inner  door. 

The  air  comes  Into  the  shelter  space  through  a  duct  over  the  inner  door 
and  is  exhausted  out  a  vent  located  at  the  rear.  The  air  intake  can  be 
closed  and  the  exhaust  vent  dropped  down  to  give  an  air  tight  closure. 
Tsinks  are  provided  for  water  and  gasoline  storage,  the  gasoline  being 
used  to  operate  an  emergency  power  unit. 

The  soil  surrounding  the  shelter  was  of  two  types  as  shown  in  Fig. 

3.  Construction  specifications  called  for  backfilling  up  to  a  level  8  ft 
from  the  floor  with  select  fill  compacted  to  a  specific  density.  The  con¬ 
tractor  chose  to  use  pea  gravel  for  this  select  fill.  The  original 

soil  removed  during  excavation  was  then  put  on  top  of  this  select  fill. 


2.2  INSTRUMENTATION 

2.2.1  Temperature  Measuring  Stations 

Upon  completion  of  the  shelter  construction,  holes  were  drilled 
through  the  walls  and  thermocouple  probes  were  Inserted  at  the  loca¬ 
tions  shown  In  Fig.  4. 

A  series  of  probes  were  located  on  the  center  line  at  U  ft,  8  ft, 
and  11  ft  from  the  floor,  at  the  top  of  the  arch  and  one  in  the  floor  at  the 
center.  An  additional  set  was  located  midway  between  the  front  and  the 
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Fig.  1  Cut-away  Diagreon  of  Shelter 
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Fip,.  2  Shelter  Entrance  Showing  Blast  Bulkhead,  Exterior 
Door,  Mushroom  Shaped  and  Intake  Vent  Covers 


k 


SURFACE  AREA  OF  SHELTER  IN  CONTACT 
WITH  THE  VARIOUS  SOILS 

SOIL  CONTACT  AREA 

SELECT  FILL  (SF)  1160  SO  FT 

BACK  FILL  (BF)  1175  m  „ 

FLOOR  (F)  1200  "  " 


Fig.  3  Soil  Conditions  Surrounding  the  ITODL  Shelter 
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B45T 


F- FRONT  SECTION  L- LEFT  SIDE 

C-CENTER  LINE  T-TOP 

B- REAR  SECTION  R- RIGHT  SIDE 

NUMBERS  DESIGNATE  THE  HEIGHT 
ABOVE  THE  FLOOR,  FOR  EXAMPLE- 
LF-B-LEFT  FRONT  8FT  FROM 
FLOOR.  RC-n- RIGHT  CENTER  11  FT 
FROM  FLOOR. 

F45T  \  THERMOCOUPLE  LOCATED 
aB45Tj  IN  CORNER  AT  45“  ANGLE. 


Fig.  L  Temperature  Measuring  Stations 
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center  line  and  another  one  midway  between  the  rear  and  the  centerline  at  the 
8  ft  level.  A  single  probe  was  located  in  the  back  wall  at  the  8  ft 
level.  The  thermocouples  at  each  station  consisted  of  five  couples, 
one  on  the  inside  of  the  metal  surface,  one  at  the  outside  svirface  and 
one  each  at  1  ft,  2  ft  and  3  ft  into  the  surrounding  fill. 

The  external  surface  and  the  thermocouples  at  1  ft,  2  ft  and  3  ft 
from  the  shell  of  the  shelter  were  positioned  on  thin-wall,  half -inch 
diameter  phenolic  resin  tubes  that  were  pushed  through  holes  in  the 
shelter  wall  (Fig.  5)*  The  thermocouple  leads  were  secured  to  the  out¬ 
side  of  the  plastic  tubes  at  proper  locations.  After  positioning  the 
four  couples  on  each  tube, one  end  was  sealed  and  the  tube  then  filled 
with  sand.  A  half  inch  pipe  with  a  loosely  fitting  steel  point  in  the 
end  was  driven  into  the  soil  at  the  desired  locations  and  the  probe 
slipped  into  place.  The  pipe  was  then  pulled  out  over  the  lead  wires, 
leaving  the  steel  point  and  the  probe  in  the  soil. 

The  internal  surface  temperature  thermocouples  were  attached  to 
the  wall  with  masking  tape.  Wet  and  dry  bulb  temperatures  were  taken 
a^  the  entrance,  center  of  shelter  and  at  the  exhaust.  The  thermo¬ 
couple  leads  were  all  taken  out  through  an  opening  in  the  top  of  the 
shelter  near  the  exhaust  vent  and  connected  to  continuous  strip  chart 
recorders.  Continuous  temperature  recordings  were  made  with  an  auto¬ 
matic  printout  of  each  temperature  every  8  minutes. 

During  the  simulated  occupancy  tests,  wet  bulb  thermocouples  were 
located  at  the  entrance  to  the  shelter,  inlet  over  the  inner  door,  cen¬ 
ter  of  the  room  and  at  the  exhaust.  The  evaporation  at  these  thermo¬ 
couples  was  accomplished  by  a  small  electric  fan. 

2.2.2  Bunan  Simulators 


The  simulators  were  designed  to  simulate  the  latent  and  sensible* 
heat  output  of  a  human.  These  human  simulators  consisted  of  a  nichrome 
wire  heating  coil  mounted  in  the  bottom  of  a  five  gallon  paint  pail. 

The  pails  were  covered  with  cotton  flannel  to  serve  as  a  wick.  Q?he 
moisture  output  of  the  simulators  was  provided  by  dripping  water  onto 
the  top  of  the  pail.  The  water  output  was  regulated  by  a  valve  at  each 
simulator  and  the  total  flow  of  water  into  the  shelter  was  controlled 
in  the  Instrument  shack.  The  amount  of  water  to  be  added  was  varied 
with  the  temperature  at  the  center  of  the  room. 

^Sensible  heat  added  to  the  shelter  is  a  direct  addition  of  heat  to  the 
enclostu-e  by  any  cneor  all  of  the  mechanisms  of  conduction,  convection, 
and  radiation.  A  gain  of  latent  heat  is  considered  to  occur  when  there 
is  an  addition  of  water  vapor  to  the  air  of  the  enclosure. 
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SHELTER 


Fig.  5  A.  Pipe  Driven  Into  Soil  P^ior  to  Inserting  Thermocouple  Probe 
B.  Thermocouple  Probe  in  Location  in  Soil  After  Removing  Pipe 


Thirty-six  simulators  vere  located  on  each  side  In  the  approximate 
position  of  the  center  bvink  as  shown  In  Fig.  6.  . 

The  remaining  28  simulators  vere  located  In  the  front  section  at 
the  tables  as  shown  in  Fig.  7* 

2.2.3  Cerbon  Monoxide  Measuring  Equipment 

The  carbon  monoxide  measuring  instrument  was  provided  by  the 
U.  S.  Forest  Service.  This  equipment >  developed  by  the  Forest  Service 
Pacific  Southwest  Forest  and  Range  Experiment  Station,  consisted  of 
surplus  detectors  which  had  previously  been  used  as  warning  device  in 
aircraft.  These  detectors  split  the  flow  of  a  gas  sample  through  two 
cells  connected  by  a  thermopile.  One  cell  contains  activated  hopcalite* 
which  catalyzes  the  oxidation  of  carbon  monoxide.  The  resulting  tei^per* 
ature  difference  between  the  two  cells  Is  a  function  of  the  carbon  mon¬ 
oxide  concentration  In  the  air  being  sampled.  QMs  temperature  differ¬ 
ence  results  In  an  emf  being  generated  In  a  thermopile  which  Is  recorded 
on  a  strip  chart  recorder.  The  unit  is  calibrated  to  give  approximately 
100  millivolts  for  a  concentration  of  one  percent  (by  volume). 


‘^Hbpcallte  -  a  mixture  of  oxides  of  copper,  cobalt,  manganese  and  silver, 
used  In  gas  masks  as  a  catalyzer  converting  carbon  monoxide  to  carbon 
dioxide . 
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CHAPTER  3 


TEST  PROCEDTIRES 


3,1  TEMPERAIURE  TESTS 

A  total  of  five  tests  vere  conducted  as  shown  in  Table  1.  During 
the  simulated  runs  No.  1  and  2,  the  water  output  of  the  simulatorB  was 
regulated  by  a  valve  at  each  station  and  the  total  flow  of  water  into 
the  shelter  was  measured  by  a  flowrator  in  the  supply  line.  The  amount 
of  water  added  was  varied  with  the  center  rocan  temperature  in  accord¬ 
ance  with  the  values  given  in  Pig.  8.  The  total  heat  into  the  shelter 
by  the  simulators  was  maintained  at  470  Btu/hr  per  man>  which  is  an 
estimate  of  the  heat  that  would  be  generated  by  a  shelter  occupant. 

This  value  is  based  on  data  from  the  Heating  and  Ventilating  Guide  of 
the  American  Society  of  Heating,  Refrigeration  and  Air-Conditioning 
Engineers  which  give  38O  Btu/hr  as  the  heat  generated  by  a  human  seated- 
at  rest,  430  standing  at  ease,  and  760,  walking  at  2  mph. 

Preliminary  heat  balance  calculations  on  runs  No.  1  and.  No.  3  indi¬ 
cated  inaccuracies  in  the  heat  transfer  coefficient  for  soil,  because 
of  the  diurnal  temperature  variations  in  the  inlet  air.  Therefore  runs 
4  and  5  vere  conducted  with  air  vents  closed.  These  tests  were  continued 
for  a  sufficient  length  of  time  to  establish  a  constant  rate  of  heat  loss 
to  the  soil. 


3.2  METHOD  OP  PREDICTING  SHELTER  THERMAL  CONDITIONS 

Heat  transfer  in  this  shelter  is  a  complex  process  and  predicting 
the  shelter  environment  under  various  climatic  conditions  is  difficult. 
Therefore  arrangements  were  made  with  the  Mechanical  Engineering  Depart¬ 
ment  of  Oklahoma  State  University  to  determine  the  feasibility  of  using 
a  system  which  they  have  developed  and  called  a  thermal  analyzer  to 
determine  shelter  thermal  behavior  under  various  specified  conditions. 
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LBS  WATER  /  HR  ADDED  TO  SHELTER  /  MAtg 


Fig.  0  Relationship  of  Water  Emitted  by  100 
PorGons  to  Ten^rature  of  Their  Environment 
(based  on  data  from  the  ASIIVE  Guide) 


Ih 


The  thermal  analyzer  applies  the  principles  of  thermal -electrical 
analogy  vhich  is  based  on  the  similarity  of  the  equations  describing 
heat  and  electrical  conduction.  The  partial  differential  equation  for 
temperature  distribution  in  a  solid  conducting  heat  can  be  compared  to 
the  partial  differential  equation  for  the  distribution  of  electrical 
potential  produced  by  a  current  flowing  through  an  electrical  network 
with  lumped  parameters.  The  rate  at  which  thermal  energy  is  stored  can 
be  compared  to  the  rate  of  flow  of  charge  into  a  capacitor.  Therefore 
a  resistance -capacitance  electrical  circuit  serves  as  an  analogue  for 
heat  flow  in  a  slab.  In  this  analogy,  electrical  current  corresponds 
to  heat  flow,  electrical  voltage  differences  to  temperature  differences, 
electrical  resistance  correapjonds  to  thermal  resistance,  and  electrical 
capacity  corresponds  to  thermal  capacity.  The  physical  system  is  lumped 
or  sub-divided  into  geometrical  sections  of  simple  shape  and  form.  The 
thermal  properties  of  each  section  are  considered  to  be  concentrated  at 
the  central  nodal  point  of  each  section  or  volume,  and  heat  is  Imagined 
to  be  conducted  between  nodal  points  throvigh  a  network  of  fictitious 
heat-conducting  rods  of  appropriate  thermal  conductance.  Details  of 
procedures  used  in  the  development  of  this  anaiysis  are  given  in  the 
Appendix. 


3.3  CARBON  MONOXIDE  MEASUREMENTS 

Air  samples  were  drawn  continuously  by  a  vacuum  pump  from  the 
front,  center,  and  rear  of  the  shelter  during  the  human  occupancy  tests 
(Run  #3)  and  put  through  the  carbon  monoxide  measuring  equipment.  Car¬ 
bon  monoxide  samples  were  also  drawn  from  these  locations  with  the 
shelter  empty  while  operating  the  gasoline  driven  motor  generator  to 
study  its  effect  on  the  shelter  environment. 


15 


CHAPTER  4 


RESUITS  AND  DISCUSSION 


4.1  TEMPERATURE  STUDIES 

Figure  9  shovs  the  shelter  temperature  during  Run  No.  1  (simulated 
occupancy).  The  relationship  of  center  shelter  temperature  to  the  air 
temperature  at  the  entrance  to  the  shelter  and  the  air  temperature  vhere 
It  enters  the  shelter  proper >  l.e.  over  the  Inner  door,  are  shovm.  The 
effective  temperature  In  the  center  of  the  shelter  Is  also  shovm.  Ef¬ 
fective  Temperature  as  described  by  ASHRAE  Is  an  empirically  determined 
Index  of  the  degree  of  warmth  perceived  on  exposure  to  different  combi¬ 
nations  of  temperature,  humidity,  and  air  movement.  The  effective  tem¬ 
perature  index  cannot  be  measured  directly,  but  is  determined  from  diy 
and  wet-bulb  temperatures  and  air  motion  observations  by  reference  to 
an  Effective  Temperature  Chart  prepared  by  the  ASHRAE.  This  was  pre¬ 
pared  from  studies  in  which  trained  subjects  compared  the  relative 
effects  of  various  air  atmospheres  by  occupying,  alternately,  two  ad¬ 
joining  air  conditioned  rooms.  The  acceptable  effective  temperature 
range  for  normally  clothed  men  and  women,  at  rest,  during  winter  and 
summer  months  is  63  to  75.  In  an  emergency  situation  this  upper  limit 
would  undoubtedly  be  raised.  At  present  an  effective  temperature  of 
85°P  Is  considered  acceptable  for  survival  shelters. 

At  the  beginning  of  the  simulated  shelter  occupancy  test  (Run  No.  l) 
the  air  entering  the  shelter  was  at  75  degrees  and  the  center  room  tem¬ 
perature  was  5  degrees  lower.  Within  4  hours  the  shelter  air  temperature 
was  at  80°  and  never  dropjred  below  75°  for  the  duration  of  the  run.  The 
shelter  temperature  reflected  the  diurnal  variations  but  was  considerably 
damped.  After  the  first  day  of  occupancy  the  entrance  air  varied  from 
94  to  50  degrees  and  the  shelter  temperature  varied  from  90  to  75  degrees. 
The  effective  temperature  varied  from  72  to  82  degrees  in  the  center  of 
the  shelter  which  is  somewhat  higher  tten  the  optimum  effective  tempera¬ 
ture. 


The  shelter  temperature  during  Run  No.  3  (Human  Occupancy)  is  shovm 
in  Pig.  10.  The  relationship  of  the  center  shelter  temperature  to  the 
entrance  and  inlet  air  temperature  are  shovm.  The  effective  temperature 
in  the  center  of  the  shelter  is  also  shown.  It  will  be  observed  that 
the  ambient  air  temperature  at  the  start  of  this  test  was  approximately 
10  degrees  lower  than  in  the  simulated  occupancy  run.  The  entrance  air 
reached  a  low  of  36  and  a  maximum  of  only  67.  The  minimum  shelter 
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temperatiire  vas  70  and  a  maximum  of  82.  Although  the  temperature  record¬ 
ings  showed  a  low  of  70  In  the  center  of  the  room  the  air  entering  the 
space  dropped  to  as  low  as  45  degrees.  Due  to  poor  inlet  air  distribu¬ 
tion  there  were  cool  drafts  in  the  front  portion  of  the  shelter.  To 
minimize  the  air  intake  problem  the  blower  was  turned  off  for  short 
periods  during  the  coolest  part  of  the  night.  This  procedure  was  fol¬ 
lowed  for  the  first  6  days.  During  the  sixth  day  the  occupants  felt 
that  they  had  corrected  the  difficiilties  by  bending  the  baffles  in  the 
inlet  duct;  however,  it  will  be  observed  on  Fig.  10  that  after  the 
sixth  day  there  was  a  coincidental  rise  in  the  inlet  air  temperature. 

The  low  entrance  air  temperature  for  the  first  6  days  was  36  degrees 
and  on  the  7th  and  8th  day  it  was  48  and  42  respectively.  After  the 
first  day  of  occupancy  the  center  room  temperature  fluctuated  from 
70°  to  82°  and  the  effective  terai^erature  varied  from  Si-®  to  74°. 

Although  Run  No.  2  was  too  short  to  give  any  significant  tempera¬ 
ture  gradients  through  the  soil  the  importance  of  the  ventilating  air 
in  removing  moistrire  from  the  shelter  was  demonstrated  Immediately  after 
the  shelter  was  closed  down.  Within  I-I/2  hrs  after  shutting  off  the 
air  supply,  molst^rre  was  condensing  on  the  ceiling  and  within  2  hrs  the 
entire  metal  surface  was  covered  with  moisture  which  started  running 
down  and  collecting  at  the  bottom.  Vfithin  10  to  15  min  after  turning 
on  the  blower  at  the  conclusion  of  the  run,  it  was  comfortable  near  the 
front  of  the  shelter  and  within  30  to  4o  min  it  was  comfortable  through¬ 
out  the  shelter. 


Analysis  of  Heat  Losses  in  Shelter 


The  heat  put  into  the  shelter  by  either  the  simulators  or  human 
occupants  must  be  expelled  through  the  shelter  walls  or  lost  to  the 
ventilating  air.  Therefore  the  total  heat  input  to  the  shelter  equals 
the  heat  lost  to  the  air  plus  the  heat  lost  to  the  soil,  or 


Q  total  = 


"‘air 


+  Q, 


soil 


(1) 


In  the  simulated  occupancy  test  the  heat  picked  up  by  the  air  consisted 
of  three  components :  q^  -  the  heat  expended  in  vaporizing  the  added 
water  and  heating  this  vapor  to  the  exhaust  temperature;  qg  -  the  heat 
expended  in  raising  the  water  vapor  in  the  inlet  air  to  the  exhaust 
temperature;  and  q3  -  the  heat  used  in  raising  the  inlet  air  to  the 
exhaust  temperature. 


\ir  ^  <13 


(2) 
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TEMPERATURE  ("F) 
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TIME  IN  hours 


Fig.  10  Temperature  Records  of  l4-Day  Human  Occupancy  of  Shelter.  Run  No 


l4-Day  Human  Occupancy  of  Shelter. 


Run  No.  3 


Shelter  human  occupancy  studies. 
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=  hg  at  tdg  -  hf  at  66°P  x  lbs  of  H2O  added 
^2  =  hg  at  td2  -  hg  at  tdj  x  x  lbs  air/hr 
^3  =  0.24  td2  -  0.24  td]^  X  lbs  dry  air/hr 

where  hg  is  Enthalpy  of  Saturated  Vapor 
ta2  Is  Exhaust  Air  Temperature 

hf  is  Enthalpy  of  Satvurated  Liquid  Supplied  to  Simulators 
tdi  is  Inlet  Air  Temperature 

S]^  is  Specific  Humidity  of  Exhaust  Air  (lbs  of  moisture/ lb  of 
dry  air) 

0.24  td2  is  the  Specific  Heat  of  the  Exhaust  Air  lbs  (dry  air/hr  = 
7730  Ibs/hr  at  I6OO  cfm) 

In  the  human  occupancy  test  qi  =  0,  since  there  is  no  water  added 
to  the  shelter  by  artificial  means. 

The  heat  lost  to  the  soil  is  expressed  as  follows: 


a  =  k  A  ^ 
^soil  ^  dx 


(3) 


where  K  is  the  Thermal  Conductivity  of  the  soil 

A  is  the  Area  of  the  Shelter  exterior  surface 

^  is  the  temperature  gradient  with  distance  from  the  shelter. 

With  the  intake  air  to  the  shelter  shut  off  during  the  calibration 
runs  Nuirtbers  4  and  5,  the  total  heat  input  is  equal  to  the  heat  lost  to 
the  soil  or 


^total  ^soil 


=  k  A  — 
^  dx 


(4) 


With  a  known  quantity  of  heat  being  put  into  the  shelter  and  the  tem¬ 
perature  gradient  in  the  soil  determined  by  measurement,  the  value  of 
k  can  be  calculated.  These  calibration  runs  (Nos.  4  and  5)  were  con¬ 
tinued  xuitll  the  heat  was  being  lost  to  the  soil  at  a  constant  rate  as 
determined  by  plotting  dt  vs  dx.  The  values  for  dt/dx  at  each  station 
after  equilibrium  was  established  are  given  in  Tables  2  and  3  for  Runs 
4  and  5  respectively.  Substitution  of  values  of  dt/dx  into  equation  4 
gave  a  value  of  1.00  for  k  on  Run  No.  4  and  1.08  on  Run  No.  5  or  an 
average  value  of  1.04. 
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The  heat  lost  to  the  air  and  to  the  soil  vas  calculated  for  each 
2-hr  inteirval  from  the  above  equations  for  both  the  simulated  and  human 
occupancy  tests.  The  results  are  given  in  Tables  4  and  5.  The  averaged 
results  for  both  runs  are  shovn  in  Table  6. 

It  will  be  noted  that  the  heat  lost  to  the  air  on  the  simulated 
run  is  68.7  i>  of  the  total  as  compared  to  69.2  5^  on  the  hxanajn  occupancy 
run.  The  calculated  total  heat  input  to  the  shelter  on  the  simulated 
run  is  only  2.7  lower  than  the  measured  input. 

The  total  average  heat  loss  during  the  human  occupancy  run  was 
43,198  Btu  per  hour  or  482  Btu  per  person.  This  result  is  compared 
with  typical  published  values  in  Table  7. 

It  is  difficult  to  extend  the  data  to  other  localities  under  vari¬ 
ous  climatic  conditions  because  of  the  large  number  of  variables  in¬ 
volved  such  as,  soil  conditions  which  affects  the  heat  loss  through 
the  walls,  diurnal  temperature  and  humidity  variations  which  affect  the 
maximum,  minimum,  and  effective  temperature  (E.T.),  and  air  flow  rate 
which  affects  the  rate  of  heat  loss  from  the  shelter  and  influences  the 
E.T.  The  data  shows  that  the  shelter  was  comfortable  for  the  entire 
l4  days  of  human  occupancy  (Run.  No.  3)  during  which  the  ambient  temp¬ 
erature  varied  between  3^  and  67°P. 

The  maximum  E.T.,  however  during  the  simulated  occupancy  test  (Run 
No.  1)  was  83  with  an  smblent  temperature  variation  between  42°  and  95°  P* 
This  E.T. is  a  few  degrees  higher  than  the  optimum  maximum  of  75* 
simulated  occupancy  test  therefore,  indicates  that  trouble  may  be  en¬ 
countered  during  the  summer,  in  areas  where  the  ambient  temperature  at 
night  does  not  drop  as  low  as  it  did  during  this  test. 


4.2  PREDICTING  THERMAL  BEHAVIOR  OP  SHEI/CER  ENVIRONMENT 

The  thermal-ele Ctrl cal  analogy  is  a  useful  tool  in  the  prediction 
of  the  thermal  behavior  of  complex  heat  transfer  systems,  but  according 
to  Oklahoma  State  University,  little  infoimation  was  available  in  the 
literature  on  its  use  in  predicting  the  thermal  behavior  of  buildings. 

The  thermal  electrics],  analogy  is  based  ution  the  identity  of  the 
fundamental  equations  of  heat  flow  within  a  rigid  body  and  tmt 
of  an  electrical  charge  in  a  non-inductive  circuit,  in  the  analysis 
made  at  Oklahoma  State  University  on  their  thermal  analyzer  the  NRDL 
shelter  interior  was  divided  into  six  air  nodes  or  cross-sections  meas¬ 
uring  8  ft  in  depth.  Node  or  section  No.  1  was  at  the  front  of  the 
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•TABUS  5  (Cont’d) 
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table;  7 

Metabolic  Rates*  of  Mea  at  Different  Activities 


Activity 

Average  Metabolic 
Rate 

SBJ  per  Hour 

Sleeping 

255 

Awake,  quiet 

300 

Seated,  at  rest 

380 

Standing,  at  ease 

430 

Normal  Shelter  Living 

482 

(Determined  In  these  studies) 

Walking,  2  mph 

760 

Walking,  k  mph 

1400 

*A11  rates  except  Konnal  Shelter  Living  are  frcm 


ASHBAE  Guide  ani  Data  Book. 
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shelter  and  No.  6  at  the  rear.  A  complete  report  of  this  analyses  is 
presented  In  the  Appendix.  Temperature  variations  of  Air  Node  No.  3 
as  determined  by  the  thermal  analyzer  is  shcrvm  in  Pig.  11  compared  to 
the  measured  temperature  in  the  center  of  the  shelter  during  the  human 
occupancy  run.  It  vill  be  noted  that  the  temperature  variations  as 
determined  by  the  thermal  analyzer  are  10  to  15  lower  than  the  measured 
temperatures.  This  was  a  preliminary  study  and  it  is  believed  that 
additional  work  would  produce  better  results.  The  analysis  procedure 
offers  the  possibility  of  an  inexpensive  and  rapid  method  of  determining 
xmderground  shelter  environment  for  a  wide  variety  of  climatic  conditions. 


4.3  CARBON  MONOXIDE  MEASUREMNTS 

During  the  htmian  occupancy  tests  the  carbon  monoxide  concentration 
was  less  than  10  ppm  from  0100  to  1000  each  day.  After  1000  each  day, 
the  concentration  increased  over  a  period  of  1  to  2  hours  to  60  to  80 
ppm  where  it  remained  until  l600  to  l800.  Thereafter,  the  concentration 
decreased  gradually  to  the  dally  low  at  about  0100.  The  increase-  in 
carbon  monoxide  concentration  can  be  explained  by  the  increased  amount 
of  smoking  during  the  daytime  hours.  The  Wavy  tolerance  for  continued 
breathing  of  carbon  monoxide  is  150  ppm.  At  no  time  during  the  occupancy 
test  was  this  level  closely  approached. 

Separate  tests  were  conducted  with  the  motor  generator  running  to 
determine  if  it  increased  the  carbon  monoxide  concentration  in  the  shel¬ 
ter.-  This  motor  generator,  which  is  located  in  the  entrance -way  oppo¬ 
site  the  ventilation  fan  Intake  is  intended  to  provide  emergency  elec¬ 
trical  power  for  ventilation  and  lights.  The  concentration  of  carbon 
monoxide  in  the  shelter  was  considerably  belcw  the  150  ppm  tolerance, 
while  running  the  motor  generator  with  the  blower  on  and  both  intake 
and  exhaust  vents  open.  Operating  the  motor  generator  however  with  the 
ventilating  air  intakes  closed^  the  carbon  monoxide  concentration  rose 
considerably  above  this  tolerance.  After  approximately  I-I/2  hrs  it 
appeared  to  be  leveling  off  at  38O  bo  400  ppm  as  shown  in  Fig.  12.  The 
carbon  monoxide  concentration  dropped  to  a  safe  value  within  a  few  min¬ 
utes  after  opening  both  the  intake  and  exhaust  air  vents  and  turning  on 
the  blower.  Although  the  motor  generator  is  not  Intended  to  be  operated 
under  these  conditions  this  test  indicates  the  hazard  that  could  be  en¬ 
countered  if  the  shelter  occupants  had  not  received  proper  instructions. 
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Fig.  11  Measured  Temperature  of  Shelter  Compared  to  Tenrperature  by  Thenno-Electrical  Analogy  for  Run  No. 


Fig-  12  Effect  of  Motor  Ctanerator  on  Shelter  Carbon  Monoxide  Concentration 
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ClIATTl®  5 


CONCLUSIONS 


The  NRDL  shelter  design  "was  proven  satisfactory  iinder  the  conditions 
of  these  tests.  The  temperature  of  the  shelter  remained  comfortable 
during  the  entire  l4  days  of  occupancy,  however  it  is  difficult  to  pre¬ 
dict  the  shelter  environment  for  other  climatic  conditions. 

The  shelter  air  temperature  reflects  the  diurnal  temperature  varia¬ 
tions  of  the  entrance  air  but  is  considerably  damped.  After  the  first 
day  of  simulated  occupancy  the  entrance  air  varied  from  50  to  94  degrees 
and  the  shelter  temperature  varied  from  74  to  90  degrees  with  an  effec¬ 
tive  temperature  variation  of  74  to  80  degrees.  During  the  human  occu¬ 
pancy  teat,  due  to  poor  inlet  air  distribution,  there  were  cold  drafts 
in  the  front  of  the  shelter  when  the  ambient  temperature  dropped  to  the 
minimum  temperature  of  36°  to  42°F.  The  maximum  entrance  air  tempera¬ 
ture  during  this  test  was  67'^.  The  air  in  the  center  of  the  shelter 
varied  from  70°  to  82°  and  the  effective  temperature  varied  from  64  to 
74  degrees. 

The  major  portion  of  the  heat  generated  in  the  shelter  was  removed 
by  the  ventilating  air.  In  both  the  human  occupancy  test  and  the  test 
run  using  hxunan  simulators,  69  ^  of  the  heat  was  carried  off  iy  the 
ventilating  air  and  only  31  $  was  lost  to  the  soil  surrounding  the 
shelter. 

The  1600  cfm  of  air  flow  through  this  shelter  reduces  the  carbon 
monoxide  from  smoking  to  a  concentration  far  below  the  Navy  tolerance 
for  submarines  of  I50  ppn,  and  it  exceeded  the  maximum  allowable  concen¬ 
tration  (mac)  established  for  industry  of  100  ppm  for  only  approximately 
1  hour.  The  industry  MAC  of  100  ppm  was  established  as  maximum  for  8 
hours  exposure.  Prom  0100  to  1000  each  day  the  concentration  averaged 
10  ppm.  The  peak  of  60  to  80  ppm  was  reached  each  day  by  1200  where  it 
remained  until  it  began  to  drop  off  between  1600  to  I8OO. 

The  location  of  the  motor  generator  will  result  in  dangerous  quanti¬ 
ties  of  carbon  monoxide  being  created  if  this  generator  is'  operated  by 
mistake  with  the  ventilating  fan  on  and  the  intake  air  vents  closed. 

Thermal-Electrical  analyzer  appears  to  offer  practical  method  of 
predicting  shelter  environment. 


CHAPTER  6 


RECO^WBNDATIONS 


It  is  reconmended  that: 

(1)  tl5ie  position  of  the  motor  generator  be  changed  or  a  safety 
measure  be  devised  to  prevent  it  from  being  operated  at  conditions 
other  than  those  for  which  it  is  intended. 

(2)  Changes  be  made  in  the  inlet  ducting  to  eliminate  cold  drafts 
at  the  center  of  the  shelter  during  cold  weather. 

(3)  Consideration  be  given  to  the  further  development  of  Thermal- 
Electrical  Analogy  as  a  tool  to  predict  shelter  environment  over  a 
wider  range  of  climatic  conditions. 
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SYMBOLS 


C 

'p 

£ 

1 

k 

m 

q 

T 

r 

V 

Hh 

-C^ 

Subscripts 


-  Capacity  (Thermal  and  Electrical) 

-  Specific  Heat  ®  Constant  Pressure,  BiO/#®F 

-  Voltage,  Volts 

-  Current  -  Amperes 

-  Thermal  Conductivity.  Btu/hr.ft,®F 

-  Mass  Flow  Rate.  )j^/hr 

~  Quantity  of  Heat  Energy  per  Unit  Time,  Btu/hr 
•  Resistance  (Thermal  and  Electrical) 

-  Density.  #m/ft^ 

-  Temperature.  ®F 

-  Time 

-  Volume,  ft^ 

-  Capacitance  (Electrical  &  Thermal) 

-  Diode 

-  Ground  Connection 

-  Resistance  (Electrical  &  Thermal) 

-  Electrical  Circuit 

-  Thermal  Circuit 
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CHAPTEa  I 


INTRODUCTION 


The  passive  R-C  electrical  analogy  for  thermal  systems 
is  a  powerful  tool  in  the  prediction  of  the  thermal  behavior 
of  complex  heat  transfer  systems,  and  many  investigations 
using  this  method  oan  be  found  in  the  literature.  However, 
very  few  appllcatlone  to  the  study  of  the  thermal  behavior 
of  buildings  Can  be  found. 

This  paper  presents  a  thermal  analysis  using  the 
electrical  analogy  to  determine  the  temperature  variation 
within  a  one  hundred  man  underground  shelter.  The  system 
studied  ia  a  Navy  Ammunition  hut  of  steel-arch  type  structure, 
48  feet  long,  with  an  arch  radius  of  12.5  feet,  which  is 
covered  by  a  minimian  of  three  feet  of  earth. 

This  structure  was  built  by  the  United  States  Naval 
Radiological  Defense  Laboratories,  (USNRDL),  approximately 
forty  miles  from  San  Francisco,  California,  where  physical 
and  analytical  studies  of  this  system  are  in  progress. 

These  studies  are  directed  toward  determination  of  environ¬ 
mental  conditions  within  an  underground  shelter,  occupied 
by  one  hundred  men  for  extended  periods  of  time.  Rough 
drafts  of  these  studies  we^’e  presented  to  the  Mechanical 
Engineering  Department  of  the  Oklahoma  State  University  by  a 
representative  of  USNRDL.  In  the  analysis  presented,  the 
system  and  boundary  conditions  are  the  same,  to  the  extent 


39 


that  it  is  possible,  as  in  the  USNRDL  tests.  This  permits 
the  results  of  this  analysis  to  be  compared  with  those  of 
the  USNRDL  tests. 

In  an  analysis  such  as  this,  the  main  problems  are  not 
associated  with  the  thermal-electrical  analogy  nor  with  the 
thermal  analyzer.  It  is  one  of  representation.  The  actual 
physical  conditions  and  boundary  conditions  must  be  represented 
by  an  equivalent  thermal  network  or  circuit.  In  accomplishing 
this,  the  physical  system  is  lumped  or  subdivided  into  a  number 
of  geometrical  subvolumes.  The  thermal  properties  of  each 
subvolume  is  considered  to  be  concentrated  at  the  central 
nodal  point  of  each  subvolume,  and  heat  is  imagined  to  be 
conducted  between  nodal  points  through  a  network  of  fictitious 
heat-conducting  rods  of  appropriate  thermal  conductance.  It 
should  be  apparent  that  the  lumping  of  the  physical  system 
and  the  choice  of  values  of  properties  assigned  to  nodal 
points  will  greatly  influence  the  accuracy  of  such  an  analysis. 
Methods  must  also  be  devised  for  the  representation  of  boundary 
conditions  and  any  other  parameters  which  might  affect  the 
thermal  behavior  of  the  system.  Chapter  II  presents  the 
equivalent  thermal  network  that  represents  the  system  and  this 
is  related  to  an  electrical  network  in  Chapter  III  by  means 
of  the  thermal-electrical  analogy.  The  electrical  circuit  is 
then  programmed  on  the  thermal  analyzer  which  leads  to  an 
analysis  which  is  interpreted  in  terms  of  thermal  properties. 
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CHAPTER  II 


THERMAL  CIRCUIT 

In  an  analysis  such  as  this,  it  is  necessary  to  idealize 
or  simplify  the  problem  so  that  it  may  be  represented  by  a 
manageable  thermal  circuit.  In  doing  this  the  following 
assumptions  were  made. 

(1)  The  shelter  is  considered  to  be  made  up  of  two  ends, 
floor,  and  aiding,  all  of  which  represent  paths  of 
heat  transfer  from  the  enclosure. 

(2)  All  thermal  energy  transferred  through  the  structure 
is  unidirectional.  This  is  to  say  that  there  is 

no  heat  transferred  between  the  paths  mentioned 
in  (1). 

(3)  Physical  conduction  paths  are  replaced  by  sections 
or  lumps  which  represent  the  thermal  properties  of 
the  path. 

(4)  The  thermal  properties  of  the  lumps  are  constant 
and  are  mean  values  over  the  temperature  range 
encountered, 

(5)  Net  radiation  exchange  within  the  enclosure  is 
negligible,  except  that  which  is  included  in  the 
conductance  from  the  men  to  the  walls, 

(6)  Solar  radiation  input  into  the  system  is 
considered  negligible. 
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(7)  Phase  changes,  such  as  condensation  on  the  walls 
are  not  considered. 

(8)  Contact  resistances  through  the  conduction  paths 
are  negligible. 

The  interior  of  the  enclosure  is  divided  into  six  nodes 
as  pictured  in  Figure  1,  The  paths  of  heat  transfer  from 
these  nodes  are  shown  in  Figures  14>  15  and  16,  in  the  Appendix. 


Figure  1.  Overall  Schematic  of  the  Enclosure 
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The  walls,  adjacent  rock,  and  earth  are  assumed  to  form 
a  series  resistance  to  the  radial  exchange  of  energy  between 
the  interior  nodes  of  the  shelter  and  its  surroundings.  The 
floor  and  adjacent  earth  also  represent  series  resistances  to 
the  energy  exchanges.  In  addition,  the  end  walls  of  the 
enclosure  represent  additional  conduction  paths  to  be  added 
to  each  end  node. 

An  interior  node  can  be  represented  by  the  network  as 
shown  in  Figure  2.  The  resistances  represent  thermal  resistance, 
and  the  capacitors  represent  thermal  capacity  of  each  lump  in 
the  conduction  path.,,  Figures  14,  15  and  16,  will  aid  in  the 
understanding  of  the  physical  picture. 


Conduction  Path  Between 
Internal  Node,  Floor  and  Earth 


Internal  Node 


Conduction  Path  Between 
Internal  Node,  Siding  and  Earth 


Figure  2,  Thermal  Circuit  Representation  of  an  Interior  Node 


The  boundary  conditions  imposed  on  the  thermal  system 
consist  of  a  heat  input  into  the  shelter  to  be  maintained 
at  470  B/hr-man  for  one  hundred  men,  and  1600  CFM  of 
ventilating  air  forced  in  at  one  end  of  the  shelter. 

The  one  hundred  men  can  be  represented  in  the  thermal 
circuit  by  a  heat  generation  or  energy  input,  “q"  into  each 
internal  node  of  the  enclosure. 

The  ventilating  air  presents  a  problem  in  that  there  is 
no  analogy  such  as  exists  for  the  conduction  of  heat. 
Physically  the  problem  is  as  follows t 
The  ventilating  air  la  assumed  to  vary  between  40^F 
and  60°F,  while  the  enclosure  temperature  is  expected  to  be 
ten  to  fifteen  degrees  above  the  mean  temperature  of  the 
ventilating  air.  The  flow  of  cooler  air  through  the 
enclosure  establishes  a  temperature  gradient  with  the 
temperature  increasing  in  the  direction  of  the  air  flow. 

This  gradient  affects  a  flow  of  energy  in  a  direction 
opposite  to  that  of  the  air  flow.  This  is  shown  schemat¬ 
ically  in  Figure  3.  Another  important  consideration  is 
that  the  temperature  of  a  particular  node  is  affected  by 
conditions  upstream  of  the  node  only. 

_ Air  Flow _ w 

O  O  O  O  O  0  0 

(0)  (1)  (2)  (3)  (4)  (5)  (6) 

_ Heat  Flow 

Figure  3*  Schematic  Showing  Direction  of  Heat  and  Air  Flow 

in  the  Enclosure 
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The  incominp,  air  impresses  its  temperature  or  potential 
on  its  entrance  point  which  will  be  designated  node  (0).  The 
increase  of  energy  of  the  air  moving  from  (O)  to  (l)  can  be 
described  by  -■  If  a  resistance  is  defined  as 

R  “  l/mCp,  then  Figure  4  would  be  an  analogous  thermal  circuit* 
The  symbol -j  |  |  represents  the  potential  of  the  incoming  air. 


Figure  4.  Ventilating  Air  Analogy  Applied  to  the  First  Air  Node 

In  considering  node  (2)  which  is  only  dependent  on  node 
(1),  while  node  (1)  is  not  dependent  on  node  (2),  we  haves 
q2  =  mCp  (T2  -  Ti) 
or 

^2  "  "‘'=p  ^"^2  -  to)  -  mCp  (Ti  -  to) 

and 

q2  =  mCp  (t2  -  to)  -  q^^ 

This  is  represented  schematically  in  Figure  5. 
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Figure  5.  Ventilating  Air  Analogy  Applied  to  First  and 

Second  Air  Nodes 

The  symbol — —  in  Figure  5  is  an  electrical  symbol 
representing  a  diode.  The  diode  used  as  shown,  prevents 
the  flow  of  energy  from  (1*)  to  (1).  This  is  a  necessary 
condition  in  that  node  (l)  is  independent  of  any  point 
down  stream  of  it.  Because  the  resistances  in  Figure  5  are 
equal,  the  potential  of  node  (!')  will  equal  that  of  node  (l), 
It  follows  that  is  dependent  only  upon  R,  t2  and  t^. 

The  circuit  is  now  evaluated  to  determine  if  it  is 
analogous  to  the  physical  conditions  it  represents. 

Physically  the  exchange  of  energy  at  node  (2)  is: 

q2  out  =  n'Cpt2 

q2  in  = 
and 

q2  net  =  mCp(t2  -  tj|^) 

It  is  concluded  that  the  physical  conditions  at  node  (2) 
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are  satisfied  by  the  thermal  circuit  of  Figure  5.  It  is 
pointed  out  again  that  the  use  of  the  diode  prevents  the 
flow  of  energy  from  (!•)  to  (1),  An  energy  balance  on 
node  (1)  in  Figure  5  would  be  identical  to  an  energy 
balance  made  on  node  (l)  in  Figure 

The  circuit  for  all  six  nodes  is  shown  in  Figure  6« 


Now  evaluating  the  enclosure  physically; 

Energy  into  enclosure  from  air  flow 
Energy  out  of  enclosure  from  air  flow  qout” 

The  net  flow  out  is  'lnet“ 

Summing  the  flow  of  energy  from  the  enclosure  using  the 
thermal  circuit; 

a  mCp  (t^  -  tj) 
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qj  =  mCp  (tj  -  t^) 

%  =  %  H  -  ^3) 

«  mCp  (t3  -  t2) 
qj  "  "iCp  (t2  -  ^1) 

qi  "  mcp  (tjL  -  tQ) 
q  =  mCp  (t^  -  to) 

Therefore  this  circuitry  ie  truly  analogoue  to  the 
physical  conditions  as  established. 

Figure  7  shows  the  thermal  circuit  for  an  interior  node 
with  boundary  conditions  included. 


US 


Figure  7.  Thermal  Circuit  Representation  of  an  Interior  Node 
with  Boundary  Conditions  Included 


CHAPTER  III 


ANAUGOUS  THERMAL  CIRCUIT 

The  thermal  electric  analogy,  which  ia  baaed  upon  the 
identity  of  the  fundamental  equationa  of  the  flow  of  heat 
within  a  rigid  body  and  that  of  a  charge  in  a  non  inductive 
circuit,  ia  well  known.^  Figure  8  ia  therefore  a  circuit 
diagram,  representing  the  electrical  as  well  aa  the  thermal 
circuit.  The  scale  factors  or  ratios  derived  to  complete 
the  analogy  are  shown  in  Table  I. 


TABLE  I 

SCALE  FACTORS  USED  IN  THE  ANALOGY 


Quantity 

Units 

Scale  Factors 

Thermal: 

Electrical: 

mmnm 

Value 

Time 

Hrs, 

Sec. 

1  hr/aec 

Capacity 

Btu/°F 

Farads 

CjICq 

10®  b/°f 

farad 

Resistance 

°F/Btu/hr 

Ohms 

10-®hr°F/B 

Potential 

Op 

Volta 

r/f 

1  op/volt 

Rate  of  Energy 
Transfer 

Btu/hr 

Coulombs/sec 

or 

Amperes 

911 

10®  Btu/hr 
amp 

^The  writer  has  written  a  paper  on  this  subject. 
Bibliography  Reference  No,  8, 
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Of  the  first  three  ratios  in  Table  I,  thermal-electrical 
analogy  permits  the  selection  of  two.  The  third  one  being 
fixed  by  the  selection  of  the  other  two.  The  relating 
equation  is  ^  ^7//?e 

It  was  desirable  in  this  analysis  to  choose  a  time  ratio 
wherein  days  would  be  sealed  down  to  a  matter  of  minutes. 

For  convenience,  one  hour  of  thermal  time  was  set  equal  to 
one  second  of  eleotrioal  time.  In  choosing  the  capacity  and 
resistance  scale  factors,  the  limitations  were  the  magnitude 
of  the  electrical  resistors  and  capacitances  available.  It 
was  also  desirable  to  keep  capacitor  leakage  at  a  minimum. 

The  choice  of  these  scale  factors  was  therefore  a  compromise. 

Of  the  last  two  ratios  in  Table  I,  the  selection  of  one 
would  determine  the  other.  The  relating  equation  for  this  is; 

’  9//  ^^7//^e 

Since  had  been  set,  the  choice  of  T/E  and  q/i  was 

again  a  compromise  between  selecting  values  of  current  of  an 
order  of  magnitude  greater  than  that  of  the  anticipated 
leakage  current  and  voltages  that  were  readily  available. 

The  boundary  conditions  imposed  upon  the  thermal  circuit 
must  be  duplicated  with  regard  to  the  electrical  circuit.  As 
mentioned  previously,  the  diurnal  temperature  variation  of  the 
ventilating  air  was  assumed  to  be  between  40®F  and  60®F. 

Having  selected  T/E  equal  to  l°F/volt,  this  was  represented 
by  a  Donner  Analog  Computer,  programmed  to  produce  an  output 
voltage  which  varied  sinusoidally  between  40  and  60  volts, 
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with  a  2k  second  period.  According  to  the  time  ratio,  24 
hours  is  represented  by  24  seconds. 

The  input  from  the  one  hundred  men  is  47j000  B/hr, 

For  each  air  node  this  amounts  to  7,833  B/hr,  or  in  the 
electrical  circuit  78.3  yt/  amps.  This  was  accomplished 
with  a  400  Volt  D.C.  power  supply  connected  by  a  large 
resistance  (approx.  4  meg._^'\_  )  to  each  air  node.  These 
resistances  were  adjusted  while  the  problem  was  being  run, 
so  that  proper  current  flow  was  achieved.  The  resistances 
labeled  2R,  31^,  4R  and  5H,  were  also  adjusted  during  running 
to  insure  that  the  desired  voltage  duplication  was  achieved. 

Figure  9  and  10  are  photographs  of  the  test  setup.  The 
problem  board  upon  which  the  electrical  circuit  is  arranged, 
is  the  panel  on  the  right  side  in  Figure  10,  The  capacitor 
leads  are  seen  on  the  left  in  the  same  photograph. 
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CHAPTER  IV 


MEASUREMENTS 


During  the  course  of  this  study,  voltage  measurements 
were  taken  at  every  node  in  the  circuit.  These  measurements 
were  recorded  on  a  Sargent  Recorder  which  received  its  input 
from  the  output  of  a  very  high  Impedence  voltage  measuring 
instrument  or  electrometer.  The  electrometers  are  an  integral 
component  of  the  Mechanical  Engineering  Department  thermal 
analyzer.  This  data  was  not  presented  in  this  report  because 
it  was  taken  only  as  a  means  to  determine  if  the  circuit  was 
performing  properly. 

Data  presented  in  this  report  is  in  the  form  of  graphs 
and  is  the  temperature  time  variation  of  the  two  middle  air 
nodes.  The  lower  left  hand  portion  of  the  curves,  (Figures 
11  and  12),  are  the  readings  at  the  middle  nodes  prior  to 
the  entrance  of  the  hundred  men  into  the  shelter.  What 
appears  to  be  a  discontinuity  in  the  curve  is  the  application 
of  the  input  representing  the  one  hundred  men.  This  was 
included  so  that  the  temperature  buildup  with  the  entrance 
of  the  men  could  be  shown.  Values  for  the  lower  left  hand 
portion  of  the  curve  should  not  be  considered  valid,  since 
the  ventilating  air  analogy  was  based  upon  the  assumption 
that  the  ventilating  air  was  of  lower  temperature  than  the 
air  inside  the  shelter. 

Measurements  of  the  leakage  resistance  of  each  node  in 
the  circuit  was  recorded.  Current  leakage  in  the  thermal 
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analyzer  represents  an  effective  loss  of  energy  from  the 
thermal  system  and  is  therefore  a  source  of  error.  Although 
this  is  undesirable^  it  must  be  realized  that  it  exists. 
Measurement  of  the  leakage  resistance  was  accomplished  by 
charging  the  capacitance  on  each  node  to  60  volts  and 
observing  the  time  for  the  node  to  drop  to  a  given  potential. 


56 


9( 


Figure  11.  Temperature  Variation  of  Air  Node  No.  4  Compared 

with  USNRDL  Data 
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Figure  12,  Temperature  Variation  of  Air  Node  No.  3  Compared 

with  USNHDL  Data 
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Figure  13.  Comparison  of  Incoming  Air  Temperature 
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CHAPTER  V 


RESULTS  AND  CONCLUSIONS 

The  results  of  this  study  are  the  curves  shown  in 
Figure  11  and  Figure  12.  Also  shown  on  these  figures  is  a 
plot  of  the  data  representing  the  center  room  temperature 
during  the  USNRDL  tests.  Comparison  of  the  center  room  data 
can  best  be  made  when  our  bovmdary  conditions  make  the  best 
possible  fit  with  actual  conditions  as  experienced  by  USNRDL. 
Figure  13  is  a  plot  of  the  actual  temperature  variation  of 
the  ventilating  air  during  the  USNRDL  test  compared  with  ^ 
our  assumed  temperature  variation. 

The  difference  in  the  results  of  this  study  and  the 
USNRDL  teat  is  of  the  order  of  to  1C®  and  is  very  good 
when  one  is  considering  a  complicated  heat  transfer  system. 

When  a  differential  equation  is  replaced  by  a  finite  difference 
equation,  errors  result.  These  errors  are  difficult  to 
determine,  but  for  the  data  taken  in  this  study  they  are 
estimated  to  be  negligible,  especially  with  respect  to 
leakage  errors. 

For  the  two  nodes  of  interest,  namely  the  two  middle 
air  nodes,  the  average  leakage  resistance  in  the  conduction 
paths  from  these  nodes  is  approximately  230  meg  ohms.  This 
corresponds  to  a  leakage  current  of  approximately  1.5  amps. 
This  is  very  small  compared  to  the  78  yO  amps  into  each  of 
these  nodes.  It  is  also  small  compared  to  the  current  which 
is  flowing  out  of  the  nodes  representing  the  ventilating  air. 
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which  is  of  the  order  of  magnitude  of  100  n  amps . 

A  more  accurate  comparison  with  the  USNRDL  tests  is  not 
possible  because  the  writer  does  not  have  detailed  information 
as  to  the  instrumentation  and  techniques  used  in  temperature 
measurements  during  the  USNBDL  tests.  The  location  of 
thermometers  and  thermocouples  with  relation  to  heat  sources 
and  proper  shielding  Is  very  important.  Any  deviation  from 
the  proper  technique  could  easily  account  for  the  difference 
shown .  -x 

In  a  last  analysis,  it  may  be  stated  that  a  thermal 
circuit  solution  can  only  be  as  good  as  the  representation 
of  the  actual  physical  system  and  the  ability  to  evaluate 
the  circuit  parameters  and  boundary  conditions.  This 
study  has  demonstrated  the  value  of  the  thermal  analyzer  in 
an  economical  and  rapid  solution  of  a  complex  heat  transfer 
system  •. 


^Author's  note:  Mr.  Othling  did  not  have  complete  experi¬ 
mental  details  at  the  time  this  report  was  prepared. 
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APPENDIX  * 


COMPUTATION  OF  THE  PARAMETERS 


Thermal  Resiatances  and  Capacities  in  the  Conduction  Path 
Throuf;h  the  Siding  (Figure  14) 

Material  8  gage  sheet  steel 

Thickness  .164A  inches 


k 


m 

27  Btu/hr.ft.°P 
.112  Btu/#°F 


Properties  of  the  earth  adjacent  to  the  shelter i 


e 

k 


108  #/£t^ 

.705  Btu/hr.ft„°F 


Cp  ,26  Btu/i5^°F 


The  earth  properties  are  those  used  in  the  USNRDL 
analytical  study  and  are  weighted  values  determined 
according  to  procedures  used  by  the  Army  Corps  of 
Engineers. 

Properties  of  the  Air  Space: 

/  13.7  ft^/# 

Cp  .24  Btu/#OF 


^Appendix  to  Othling's  original  document. 


For  reflective  surfaces  with  an  effective  emissivity  of 
0.05,  the  "Heatingj  Ventilating  and  Air  Conditioning  Guide" 
recommends  an  inside  surface  film  coefficient  of  0.8. 

This  includes  convection  and  radiation  effects  and  is 
independent  of  the  position  of  the  surface  and  the 
temperature  difference.  Even  though  this  particular 
choice  of  "h"  is  subject  to  discussion,  it  is  of  minor 
importance,  in  that  it  represents  a  small  percentage  of 
the  total  resistance  to  heat  flow  through  the  walls  of 
the  shelter. 
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h  =  .8  Btu/hr.ft.°F  is  the  value  used  in  the  USNRDL 
tests  and  likewise  will  be  used  in  this  study. 

o  =  _Z_  r  _ _ 

/)/!  .S 

=  3.98x10“^  hr°F/Btu 

r  r  jOCfi/  - 

r  ^  f3.7xi 


=  34.A  Btu/°F 


”2  g  /^Aj.  ^~A'A  ' ' 

=  =  .OT*ioArOr/Btu 


=  236  Btu/°F 


«3 


=  .00434  hr®P/Btu 


C 


3 


-  /as  X.UX^/^^'K*)i9 
'S 


a  1896  Btu/°P 


Thermal  Resistances  and  Capacities  in  the  Conduction  Path 
Through  the  Ends  of  the  Shelter  (Figure  15) 

Material  7  gage  sheet  steel 

Thickness  .1793  inches 
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Figure  15.  Conduction  Path  through  End  Walls 
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k  .75  Btu/hr.ft/P 

Cp  .2  Etn/f? 


nnn 

)Air  Node 

t«7 

> 

1  I 

tRg 

^  Floor 

1 

< 

■ 

ilH^ 

.£ _  e. 

Earth 

3 

i 

- 

■■i 

Figure  16.  Conduction  Path  through  Floor 


R." 

%  --  -it 

Cg  = 

"9"  iit 

C,  = 


-  1  =  .00625  hr°F/Btu 

7Sc25x§ 

=  1.5  =  .833xlO~3hr°F/Btu 

12x.75x25x^ 

=400x. 2x2 5x8x3  =  4,000  Btu/°F 
12 


~  1  =  .00708  hrOF/Bbu 

7705x2 5x^ 

=108x. 26x2x2 5x8  =  11,200  Btu/°F 


Resistance  R  =  l/jnCp  defined  by  ventilating  air  analogy; 
Flow  rate  =  1600  C.F.M. 

R  =  1/mc  =  13.7  =  .594xlO"3hr®F/Btu 

P  i60Qx. 24x60 
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TABLE  II 


SUMMARY  OF  RESISTANCES 


Symbol 

Re 

399x10“^ 

39.9/’.y-\_  Convective  &  Side  Walls 

.00434 

.434  meg^^  Side  Walls  &  Earth 

'4 ’'"5 

.0509 

5.09  meg^-N.  Convective  &  End  Walls 

B5  i'  R6 

.00578 

.578  meg../~\.  End  Walls  &  Earth 

By  l'  Bg 

.00708 

.708  meg.r^  Convective  &  Floor 

Rg  /  R9 

.00791 

.791  megy-N- Floor  &  Barth 

R  =  1/mCp 

594x10-3 

59.4 Ventilating  Air 

Xr=  Rt  “ 

10"®  hr^F/B 

% 

ohm 
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TABLE  III 


SUMMARY  OF  CAPACITANCES 


3yinbol 

Ce 

Cl 

34.4 

.34^/ 

Internal  Air 

236 

2.3^ 

Side  Wall 

<=3 

1896 

18.9^ 

Earth  AdJ.  to 
Side  Wall 

°5 

200 

Zyt// 

End  Walls 

1.38x10^ 

13|>/ 

Earth  AdJ.  to 
End  Wall 

Ca 

4,000 

kOyi// 

Floor 

<=9 

11,200 

115/// 

Earth  AdJ.  to 
Floor 

=  °t  - 

10®  B/OP 

Farad 
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